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Abstract. In this paper we study fc-noncrossing RNA structures with arc-length > 3, i.e. RNA 
molecules in which for any i, the nucleotides labeled % and i + j (j = 1, 2) cannot form a bond 
and in which there are at most k — 1 mutually crossing arcs. Let S^^n) denote their number. 
Based on a novel functional equation for the generating function 5Z n>0 Sk,3( n ) z ™ > we derive for 
arbitrary fc > 3 exponential growth factors and for fc = 3 the subexponential factor. Our main 
result is the derivation of the formula S 3 3 (n) ~ 6.11170-4! 4,54920". 



1. Introduction 

RNA is believed to be central for the understanding of evolution. It acts as genotypic legislative in 
form of viruses and viroids and as phenotypic executive in form of ribozymes, capable of catalytic 
activity, cleaving other RNA molecules. This dualism gives rise to the hypothesis that RNA may 
have preceded DNA and proteins, therefore playing a key role in prebiotic evolution. In light 
of growing support of an RNA world JTS] and even RNA-based metabolisms or the prospect of 
self-replicating RNA [T5] it is the phenotypic aspect of RNA that still lacks deeper understanding. 
Despite the fact that pseudoknot RNA structures are known to be of central importance [1] , little 
is known from a theoretical point of view, for instance only recently their generating function has 
been obtained [8]. 

Let us provide first some background on RNA sequences and structures which allows us to put our 
results into context. The primary sequence of an RNA molecule is the sequence of nucleotides A, 
G, U and C together with the Watson-Crick (A-U,G-C) and (U-G) base pairing rules specifying 
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the pairs of nucleotides can potentially form bonds. Single stranded RNA molecules form helical 
structures whose bonds satisfy the above base pairing rules and which, in many cases, determine or 
are even tantamount to their function. Three decades ago Waterman et.al. pioneered the concept 
of RNA secondary structures [HI [2TJ [T71 [7] , a concept being best understood when considering 
structure as a diagram, drawing the primary sequence of nucleotides horizontally and ignoring all 
chemical bonds of its backbone. Then one draws all bonds, i.e. nucleotide interactions satisfying 
the Watson-Crick base pairing rules (and G-U pairs) as arcs in the upper halfplanc, thereby 
identifying structure with the set of all arcs. Then secondary structures have no two arcs 
(*2, where i\ < j\ and i-i < with the property i\ < 12 < ji < ji and all arcs have at least 
length 2. While the concept of secondary structure is of fundamental importance it is well-known 
that there exist additional types of nucleotide interactions [1] . These bonds are called pseudoknots 
[2"3"j and occur in functional RNA (RNAscP [H]), ribosomal RNA [TU] and are conserved in the 
catalytic core of group I introns. In plant viral RNAs pseudoknots mimic tRNA structure and in in 
vitro RNA evolution |20j experiments have produced families of RNA structures with pscudoknot 
motifs, when binding HIV-1 reverse transcriptase. 

Leaving the paradigm of RNA secondary structures, i.e. studying RNA structures with crossing 
bonds, poses challenging problems for computational biology. Prediction algorithms for RNA 
pseudoknot structures are much harder to derive since there exists no a priori tree-structure and 
the subadditivity of local solutions is not guaranteed. In fact pseudoknot RNA structures are 
genuinely non-inductive and seem to be best described by the mathematical language of vacillating 
tableaux [5J [5] • One approach of categorizing RNA pscudoknot structures consists in considering 
the maximal size of sets of mutually crossing bonds, leading to the notion of fc-noncrossing RNA 
structures [8J. This concept is very intuitive, a fc-noncrossing RNA structure has at most k — 1 
mutually crossing arcs and a minimum bond-length of 2, i.e. for any i, the nucleotides i and i + 1 
cannot form a bond. In this paper we will consider fc-noncrossing RNA structures with arc-length 
> 3. Their analysis is based on the generating function, which has coefficients that are alternating 
sums. This fact makes even the computation of the exponential growth factor a difficult task. To 
make things worse, in case of arc-length > 3 there exists no explicit formula for the coefficients, 
which could only be computed via a recursion formula. 



1.1. Organization and main results. Let Sk,3(n) (fc > 3) denote the set of fc-noncrossing 
RNA structures with arc-length > 3 and let Sfc^n) = ^^(n)!- In Section [2] we provide the 
necessary background on fc-noncrossing RNA structures with arc-length > 3 and their generating 
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FIGURE 1 . fc-noncrossing RNA structures, (a) secondary structure (with isolated labels 
3, 7, 8, 10), (b) non-planar 3-noncrossing structure 



function Yln>o ^k,3(n)z n . In Section [3] we compute the exponential factor of Sfc. 3 (n). To make 
it easily accessible to a broad readership we give an elementary proof based on real analysis and 
transformations of the generating function. Central to our proof is the functional identity of 
Lemma [5] and its generalized version in Lemma El In Section [5] we present the asymptotic analysis 
of §3,3(n), using Flajolet et.aVs singular expansions and transfer theorems [U El El HH EH ■ This 
analysis is similar to [8] but involves solving a quartic instead of a quadratic polynomial in order 
to localize the singularities. The main result of the paper is 

Theorem. The number of 3-noncrossing RNA structures with arc length > 3 is asymptotically 
given by 



6 11170 • 4' 

1-1) S 3>3 n ~ -, — -. -4.54920", 

n(n — 1) ... (n — A) 

^ere s^{n) = n( ^n% = 146 - 6807 " 3 + + O^)] ■ 

In the table below we display the quality of our approximation by listing the subexponcntial factors, 
i.e. we compare for k = 3 the quantities Sfe i 3(n)/(4. 54920)™ obtained from the generating function 
(Theorem [1]), which are the exact values and the asymptotic expressions s 3 3 (n), respectively. 
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The sub exponential factor 



n 


S 3 , 3 (n)/(4.54920)" 




n 


S 3 , 3 (n)/(4.54920)" 


s 3 , 3 W 


10 


3.016 x 10~ 4 


4.851 x 10" 


-3 


60 


3.457 x 10~ 7 


2.238 x 10" 


7 


20 


2.017 x 10~ 5 


7.884 x 10" 


-5 


70 


1.476 x 10~ 7 


1.010 x 10" 


7 


30 


3.513 x 10~ 6 


8.577 x 10" 


-6 


80 


3.783 x 10^ 8 


5.085 x 10" 


-8 


40 


9.646 x 10~ 7 


1.858 x 10" 


-6 


90 


2.154 x 10~ 8 


2.781 x 10" 


-8 


50 


5.627 x 10" 7 


5.769 x 10" 


-7 


100 


1.299 x 10" 8 


1.624 x 10" 


-8 



2. fc-NONCROSSING RNA STRUCTURES WITH ARC-LENGTH > 3 

Suppose we are given the primary RNA sequence 

ACUCAGUUAGAAUAGCCGGUC . 

We then identify an RNA structure with the set of all bonds different from the backbone-bonds of 
its primary sequence, i.e. the arcs + 1) for 1 < i < n — 1. Accordingly an RNA structure is a 
combinatorial graph over the labels of the nucleotides of the primary sequence. These graphs can 
be represented in several ways. In Figure [5] we represent a structure with loop-loop interactions in 
two ways. In the following we will consider structures as diagram representations of digraphs. A 

C A 




AC UCAGUUAGAAUAGCC GGUC 

FIGURE 2 . A 3-noncrossing RNA structure with arc-length > 3, as a planar graph (top) 
and as a diagram (bottom) 

digraph D n is a pair of sets Vd„ , Eo n , where Vd„ = {1, ■ ■ ■ , n} and Eo n C {(£, j) \ 1 < i < j < n}. 
Vd„ and Eo n are called vertex and arc set, respectively. A fc-noncrossing digraph is a digraph in 
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which all vertices have degree < 1 and which does not contain a fc-set of arcs that are mutually 
intersecting, or equivalently 

(2.1) /3 (*rij.7ri)i (*r 2 ? Jr 2 )i •• • j (*r fc ? jr k )i Vi < Va < ' ' ' < V fc < Jn ^> jV 2 < ' ' ' < • 

We will represent digraphs as a diagrams (Figure [2]) by representing the vertices as integers on a 
line and connecting any two adjacent vertices by an arc in the upper-half plane. The direction of 
the arcs is implicit in the linear ordering of the vertices and accordingly omitted. 

Definition 1. An fc-noncrossing RNA structure with arc-length > 3 is a digraph in which all 
vertices have degree < 1, having at most a k — 1-set of mutually intersecting arcs without arcs of 
length < 2, i.e. arcs of the form + 1) and + 2), respectively. Let Sk,3(n) and 5k,3(n,£) 
be the numbers of fc-noncrossing RNA structures with arc-length > 3 and those with £ isolated 
vertices, respectively. 



Let fk{n,£) denote the number of fc-noncrossing digraphs with I isolated points. We have shown 
in [8] that 

(2.2) fk(n t £)= (^jf k (n-£,0) 

x n 

(2.3) dct[V,(2:r) - I i+j (2x)]\^i 1 = £ f k (n, 0) — 

n>l 

£ n ( " ~) n 

(2.4) e^det[I i . j (2x)-I i+j (2x)}\^ 1 = (£ |-)(£ /*M)2L) = £ £/ fe (n,£) ^ . 

£>0 ' n>l ' n>l (.1=0 ) 

where I r (2x) = J2j>a j\(r+j)\ 1& * nc hyperbolic Besscl function of the first kind of order r. In 
particular we have for k = 3 



(2.5) f 3 (n,£) 



2 + z 2 —5 HI 



where C m = m * 1 ( 2 ™) is the mth Catalan number. The derivation of the generating function 
of fc-noncrossing RNA structures, given in Theorem [1] below uses advanced methods and novel 
constructions of enumerative combinatorics due to Chen et.al. and Stanley's mapping between 
matchings and oscillating tableaux i.e. families of Young diagrams in which any two consecutive 
shapes differ by exactly one square. The enumeration is obtained using the reflection principle due 
to Gessel and Zcilbcrgcr [6] and Lindstrom [llj in accord with an inclusion-exclusion argument in 
order to eliminate the arcs of length < 2 [8] . 
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Theorem 1. Let k £ N, k > 2. Then the numbers of k-noncrossing RNA structures Sk y s{n,(.) and 
Sfc,3(n) are given by 

(2.6) S k , 3 (n,£) = ]T(-l) b A(n,6)/ fc (n-26),£) 

6>0 

|n/2j r n -2b -\ 

(2.7) S M (n) = ]T(-l) b A(n,&)<^/ fe (n-2M) • 

6=0 I £=0 J 

where A(n, 6) satisfies the recursion 

(2.8) A(n, b) = A(n -1,6) + A(n - 2, 6 - 1) + A(n - 3, 6 - 1) + A(n - 4, 6 - 2) 
anrf i/ie initial conditions for eq. (|2.8|) are A(n, 0) = 1 and A(n, 1) = 2n — 3. 



3. The exponential factor 



In this section we obtain the exponential growth factor of the coefficients Sk,3(n). Let us begin by 
considering the generating function J2 n >o Sk,3(n)x n as a power series over R. Since X)n>o Sfe(n)x TI 
has monotonously increasing coefficients lim,,-^ Sfc^n)™ exists and determines via Hadamard's 
formula its radius of convergence. Due to the inclusion-exclusion form of the terms Sk,3{n), it 
is not obvious however, how to compute this radius. Our strategy consists in first showing that 
Sfc,3( n ) i s closely related to /fc(2n,0) via a functional relation of generating functions. 

Lemma 1. Let x be an indeterminate over M. and T^in) be the number of k-noncrossing partial 
matchings over [n], i.e. Tk(n) = X)m<[-J {2m) f k fi m > 0)- Let furthermore pk denote the radius of 
convergence of^2 n>0 Tk(n)x n . Then we have 

(3.1) VM< Pfe ; ^T fc (n)x" = T ^- ^/ fe (2n,0) 

n>0 n>0 




Proof. Wc will relate the ordinary generating functions of Ife(n) 

M £ T *( n )£f =£ £ (^V,(2-,0)^ = e ^det[V J (2x)-/ l+J (2x)]^ 1 

n>0 ' ra>0m<§ ^ ' 

(see eq. (|2.4p ) and /fc(2n, 0) via Laplace transforms as follows: J2 n>0 Tk(n)x n = J2 n>Q Tk(n)'^n\ 
is convergent for any whence we derive, using the Laplace transformation and interchanging 
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integration and summation 

E r *(»)*"=E r *(»)|r r r ewt'" 

n>0 n>0 J0 170 n>0 

We next interpret the rhs via eq. (|3.2|) and obtain 

/>oo 

Y / T k {n)x n = / e a;t det[/ i _ J (2a ; t)-/ i+J (2xt)] l fe -i 1 e- t di 
Interchanging integration and summation yields 

(xty 



rOO 

n>0 J ° n>0 ^ '' 



1 Z" 00 

SA(2»,0)^.(2«)K T ^)-. T ^ 

n>0 V ; 

J-^/ fc (2n,0)( T ^-) 
1 — x ^-^ 1 — x 



dt 

• I - x ' 

1 



.r ' 



n>0 

and the proof of the lemma is complete. □ 

Lemma 2. Let x be an indeterminante over K and Sfc^n) &e </ie number of k-noncrossing RNA 
structures with arc-length > 3. JTien we Ziaue the functional equation 

(3.3) £ s M( n) = r^r^r^ E o) ( r 

n>0 n>0 v 



Proof. According to Theorem [T] we have 



x — x 



n-2b 



Sfc,s(») = E i-^K^b) E fk{n-2b,i) 

6<LfJ ^=0 

= £(-l)>AM)f; f^)/*(m-26,0) 



b<[fj m=2b 
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Our goal is now to relate J2n>o ^>k,3(n)x n to the terms Tfe(n). For this purpose we derive 



^S ti3 W/ = ^^(-lfA(^) Y, 



n>0 



n>0 2b<n 



m=2b 



- 26 

i -2b 



f k (m-2b,0)x r ' 



^(-1)V X(n,b)T k (n-2b)x 

b>0 n>2b 

^(-1) V b X ( n + 2b ' b)T k (n) x' 



i-26 



b>0 



n>0 



Interchanging the summations w.r.t. b and n we consequently arrive at 



(3.4) 



£>,3(n)* n = E 



n>0 



n>0 



^(-l)V fc A(n + 26,6) 



fa>0 



T k {n)x n 



We set ip n (x) = J2b>o M n + 26, 6)a; b . According to Theorem Q] we have the recursion formula 

A(n + 26, 6) = A(n + 26 - 1, b) + X(n + 26 - 2, 6 - 1) + A(n + 26 - 3, 6 - 1) + A(n + 26 - 4, 6 - 2) 

Multiplying with x b and taking the summation over all 6 ranging from to [n/2\ implies the 
following functional equation for ip n (x), n = 1,2... 



(3.5) 



lfi n {x) = lf„-l(x) + X ■ lfi n (x) + X ■ (p„-l(x) + X 2 lf n (x) 



Eq. (|3.5|) is equivalent to 



(3.6) 



1 + X 



(p n -i(x) 1-x-x 2 



and (p (x) = V A(26, 6)x b = - 



x — x" 



since A/, = A(26, 6) satisfies the recursion formula Ah = Ae,_i + \b-2 and the initial condition 
Ao = Ai = 1. We can conclude from this that A(26, 6) is 6-th Fibonacci number. As a result we 
obtain the formula 



(3.7) 



ip„{x) = tp (x) 



l + x 



1 — x — x* 



1 



1 — X — x 2 \ 1 — x — X 



l+x 
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Substituting eq. (|3.7p into eq. (|3.4[) we can compute 



n>0 



^ l + x 2 -x 4 (l + x 2 -x 4 ) Tk{n)xT ' 

n>0 

i + j-^ rt(n) (r 



+ x 2 — x 4 



z — ' 1 — x ^— ' 1 — X 



n>0 

Via Lemma [T] we have the following interpretation of X) n >o ^fe( n ) ^ 

^2n 

jfc^/t) u 

1 — X 

n>0 ~~ n>0 

Therefore we obtain setting; x l = = .^T/ — T , 1-^—7 = ■= — 1 ^ ~^ — T and -^—r = ^ >> — T 

l-\-x z — 1 1 — x' l—x+x^+x^—x^ 1—x' 1 — x-\-x*-\-x^ — X* 

n>0 n>0 

whence the lemma. □ 



Using complex analysis we can extend Lemma [2] to arbitrary \z\ < pk, where z £ C 

Lemma 3. Let k > 2 be an integer, then we have for arbitrary z £ C with the property \z\ < pk 
the equality 

1 / . 3 \ 2n 

0.8) e s m (») - n = 1 _ z + Z 2 + z3 , Z 4 E /^( 2 "' °) r 



n>0 ti>0 



— 2 + Z z + Z d — 



In 

are 



Proof. The power series J2n>a S kA n ) zn and TZi+i^jrp- E„>o fk{2n 7 0) ( \_ Z 44'I 3 - 
analytic in a disc of radius < e < and according to Lemma [2] coincide on the interval ] — e, e[ 
Therefore both functions are analytic and equal on the sequence (^) n eN which converges to and 
standard results of complex analysis (zeros of nontrivial analytic functions are isolated) imply that 
eq. (|3.8[) holds for any z £ C with \z\ < pk, whence the lemma. □ 



Lemma [3] is the key to compute the exponential growth rates for any k > 2. In its proof we recruit 
the Theorem of Pringshcim |19j which asserts that a power series X)«>o a n z ™ with a„ > has 
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its radius of convergence as dominant (but not necessarily unique) singularity. In particular there 
exists a dominant real valued singularity. 



Theorem 2. Let k > 3 be a positive integer and r^ be the radius of convergence of the power series 

i>0 • 



£ n>0 / fe (2n,0)z 2 ™ and 



(3.9) * :[0> >S^[0, 5 -A ~ *(!-*)(! + *) 



2 J L ' 4 J ' _(*2_ I)2 + z(z 2 - + 5 ■ 

Tften i/ie power series ^2 n>0 Sk t 3{n)z" has the real valued, dominant singularity pk, which is the 
unique real solution of i9(x) = r^ and for the number of k-noncrossing RNA structures with arc- 
length > 3 holds 

(3.10) S M (n) ~ (— 

\Pk 

In Section 2] we will in particular prove that p$ « 0.21982. 



Proof. Suppose we are given r^, then r^ < ^ (this follows immediately from C n ~ 2 2 ™ via Stirling's 
formula) . The functional identity of Lemma [5] allows us to derive the radius of convergence of 
Y^ n >o Sk(n)z n . According to Lemma[3]we have 

1 / 3 \ 2n 

(3.1D £ s M( „) *- = r37T ^ T ^ e o) ( 1 _ z ; z ;; z3 , z4 ) 

n>0 n>0 

/fe(2n, 0) is monotone, whence the limit lim^oo /fe(2n, 0) A: exists and applying Hadamard's for- 
mula: lim n ^oo fk(2n,0)^n = For z e 1, we proceed by computing the roots of 



''A- 



1 — z + z 2 + z 3 — z 4 

which for rk < \ has the minimal root pk- We next show that pk is indeed the radius of convergence 
of X) n >o Sk(n)z n . For this purpose we observe that the map 
(3.12) ~ 

V2 rn 5-^2 z(l-z)(l + z) 

*:[0,-2-] — [0,-^], — _ (z2 _i )2 + z(z2 _i)_ § + r *0*)=r* 

is bijective, continuous and strictly decreasing. Continuity and strict monotonicity of •d guarantee 
in view of eq. (|3.12[) that p/., is indeed the radius of convergence of the power scries X) n >o ^k(ji)z n . 
In order to show that pk is a dominant singularity we consider £n>o Sfc,3(^)z™ as a power series 
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FIGURE 3. We display 4 poles(the corresponding 4 peaks) in 1 _ z+ z z 'i+ z >j_,4, over C. The 
picture illustrates that i? is a bijection on the interval [0, -^], which allows us to obtain 
the dominant singularity p^. 



over C. Since S^^n) > 0, the theorem of Pringsheim [19] guarantees that p k itself is a singularity. 
By construction p k has minimal absolute value and is accordingly dominant. Since S^^n) is 
monotone lim„_, 00 Sk,3{n)i exists and we obtain using Hadamard's formula 

(3.13) lim Sfc3(n)~ = — , or equivalently S/c3(n)~ [ — 

n— oo p k \p k 

from which eq. (|3.10|) follows and the proof of the theorem is complete. □ 




4. Asymptotic Analysis 



In this section we provide the asymptotic number of 3-noncrossing RNA structures with arc- 
length > 3. In the course of our analysis we derive the analytic continuation of the power series 
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Sn>o ^3,3( n ) z ™- The analysis will in particular provide independent proof of the exponential fac- 
tor computed in Theorem [2j The derivation of the subexponcntial factors is based on singular 
expansions j¥j in combination with transfer theorems. The key ingredient for the coefficient ex- 
traction is the Hankcl contours, see Figure [U Let us begin by specifying a suitable domain for our 
Hankcl contours tailored for Theorem [5J 

Definition 2. Given two numbers <fi, R, where R > 1 and < <j> < % and p £ M. the open domain 
A p (0, R) is defined as 



A domain is a A p -domain if it is of the form A p (</>, R) for some R and <p. A function is A p -analytic 
if it is analytic in some A p -domain. 



FIGURE 4. Ai-domain enclosing a Hankel contour. We assume z = 1 to be the unique 
dominant singularity. The coefficients are obtained via Cauchy's integral formula and the 
integral path is decomposed in 4 segments. Segment 1 becomes asymptotically irrelevant 
since by construction the function involved is bounded on this segment. Relevant are 
the rectilinear segments 2 and 4 and the inner circle 3. The only contributions to the 
contour integral are being made here, which shows why the singular expansion allows to 
approximate the coefficients so well. 





y 




Since the Taylor coefficients have the property 



(4.2) 



V 7 GC\0; \z n ]f{z)= 1 n [z n ]f(~) , 
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we can, w.l.o.g. reduce our analysis to the case where 1 is the dominant singularity. We use 
U(a, r) = {z G C\\z — a\ < r} to denote the open neighborhood of a in C. 

We use the notation 

(4.3) (f{ z ) = O (g(z)) as z — > p) (f(z)/g(z) is bounded as z — > p) 

and if we write /(z) = 0(g(z)) it is implicitly assumed that z tends to a (unique) singularity. 
[z n ] f(z) denotes the coefficient of z n in the power series expansion of f{z) around 0. 

Theorem 3. [3] Let a be an arbitrary complex number in C\Z<o, f(z) be a Ai-analytic function. 
Suppose r G Z>o, and f(z) = 0((1 — z) r \n(j^)) in the intersection of a neighborhood of 1 and 
the Ai-domain, then we have 

(4.4) [z n ]f(z)~K(-l) r -. -p— for some K>0. 

n{n — 1) . . . (n — r) 



We are now prepared to compute an explicit formula for the numbers of 3-noncrossing RNA 
structures with arc-length > 3. 

Theorem 4. The number of i-noncrossing RNA structures with arc length > 3 is asymptotically 
given by 

S »<»> - , "" 7 °; 4! S 4.54920V 
n(n — 1) ... (n — 4) 



Proof. Claim 1. The dominant singularity p^ of the power series X) n >o ^3, 3(^)2™ is unique. 
In order to prove Claim I we use Lemma [31 according to which the analytic function S 3 (z) is the 
analytic continuation of the power series X) n >o ^3( n ) z ™- We proceed by showing that S 3 (z) has 
exactly 12 singularities in C and the dominant singularity is unique. The first four singularities are 
the roots of the quartic polynomial P(z) = 1 — z + z 2 + z 3 — z 4 . Next we observe in analogy to our 
proof in [9] that the power series J2n>o f^i^ 71 ' 0)j/ n has the analytic continuation Hf(y) (obtained 
by MAPLE sumtools) given by 

(4.5) m = -i — g> T 16 ^ , 

Id j/2 

where P™(x) denotes the Lcgcndre Polynomial of the first kind with the parameters v = | and 



m 



T. ^>(y) has one dominant singularity at y — which in view of $(z 
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induces exactly 8 singularities of S3 (2) 



l+z 2 -z i -z+z' 



Indeed, ^(y 2 



l+z'-'-z^-z+z 3 

has the two singularities C: j3\ = \ and /?2 = — \ which produces for S 3 (z) (solving the quartic 
equation) the 8 singularities p 3 w 0.21982, C2 ~ 5.00829, Cs « -1.07392, C4 ~ 0.84581, ( 5 w 
-0.53243 + 0.1195H, Ce ~ -0.53243 - 0.11951i, C? ~ 1.10477 and ( 8 w -3.03992. The above 
values have error terms (for details for solving the general quartic equation see Section [5]) of the 
order 10 -5 , which allows us to conclude that the dominant singularity p 3 is unique and Claim 1 
follows. 



Claim 2. [9] ^(z) is (<fi, i?)-analytic and has the singular expansion (1 — 16z) 4 ln Ijzj 



(4.6) Vze A A (0,i?)nf7(^,e); = O ( (1 - 16z) 4 In 



1 - 16z, 

First A_i_ (0, i?)-analyticity of the function ^>(z) is obvious. We proceed by proving that (1 — 

16z) 4 ln ( 1 _\ 6z J is its singular expansion in the intersection of a neighborhood of and the A- 
domain Aj_ (</>, R). Using the notation of falling factorials (n — 1)4 = (n — l)(n — 2)(n — 3)(n — 4) 

16 

we observe 



f 3 (2n,0) = C n+2 C n -C* 



1 



12(ra- l) 4 (2n+ 1) 



2n 

(n - 1)4 (?i + 3)(n + l) 2 (?i + 2) 2 V n 
With this expression for f 3 (2n, 0) we arrive at the formal identity 



]ri6-"./3(2n,0); 



rs>5 



n>5 

+E 



16 



1 



12(n- l) 4 (2n+ 1) 



2n 



(n-l) 4 (n + 3)(n+l) 2 (n + 2) 2 V.™ 
4! 11 

2") , 



4! 



1 1 



(n — 1)4 7r n 



(n— 1)4 7r n 



where f(z) = 0(g(z)) denotes that the limit f(z)/g(z) is bounded for z — > 1, cq. (14.3)) . It is clear 
that 



n>5 



lim( > 

E 



16"" 



1 



12(n - l) 4 (2n + 1) 



(n- 1)4 (n + 3)(n+ l) 2 (n + 2) 2 



n>5 



16" 



1 



12(n- l) 4 (2n + 1) 



(n- 1)4 (n + 3)(n+ l) 2 (n + 2) 2 



2// 
n 

2 



4! 



1 1 



{n — 1)4 7r n 



4! 



1 1 



(n — 1)4 7r n 



< K 



for some k < 0.0784. Therefore we can conclude 



(4.7) 



^16-"/ 3 (2n,0)z" = O(^ 



4! 



n>5 



' in — 1)4 7r n 
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We proceed by interpreting the power series on the rhs, observing 



(4.8) 



Vn>5; [z 11 ] (1 - z) 4 In 



4! 



1 — 2/ (n — 1) . . . (n — 4) n 



whence ^(1 — z) 4 In y 1 -^ is the unique analytic continuation of En>5 (injlin 2 "' Using the 
scaling property of Taylor coefficients [z n ]f(z) = J n [z n ]f(— ) we obtain 



(4.9) 



'zeAj_ (<f>, R) n f/(^,e); = O ((1 - 16z) 4 In (tZ^)) 



Therefore we have proved that (1 — 16z) 4 ln( 1 _ 1 16z ) is the singular expansion of ^(z) at z — -j^, 

whence Claim 2. Our last step consists in verifying that the type of the singularity does not change 

when passing from («) to S 3 (z) = x _ z+z i+ zi _ z , *(( ) 2 )- 

Ctem 3. For z € A P3 (0, i?) n U(p 3 , e) we have E 3 {z) = o((l- ^) 4 ln( T ^)). 

To prove the claim set u(z) = 1 — z + z 2 + z 3 — z A . We first observe that Claim 2 and Lemma[3] 

imply 



S 3 (z) = O 



u(z) 



i-mC-j^-) 2 ) In 

u(z) 



1 _ 16(^1) 2 



The Taylor expansion of g(z) = 1 — 16( z u ^ ) 2 at p 3 is given by q(z) = a(p 3 — z) + 0(z — p 3 ) 2 and 
setting a « —1.15861 we compute 



1 



u(z) 



q(z) 4 hi^- 



(q(p 3 ~ j + Q(z - p 3 ) 2 ) 4 In a(p3 - z) +o( z - P3 )* 
0.83679 - 0.45789(z - p 3 ) + 0((z - p 3 ) 2 ) 

([a + 0(z - P3 )]{ P3 - z) 4 In [a+0{z J P3)]{p3 _ z) 



0(z - P3 ) 



= 0({p 3 -zf\n^— ) , 

P3 - Z 

whence Claim 3. Now we are in the position to employ Theorem [31 and obtain for S 3i3 (n) 



S 3 , 3 (n)~K' [z n ] (p 3 -z) 4 ln 



K 1 ■ 



4! 



p 3 - z J n(n — 1) ... (n — 4) \p 3 



Theorem [T] allows us to compute K' = 6.11170 and the proof of the Theorem is complete. □ 
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5. Appendix 



Let us first introduce some basic definitions for quartic equation (quartic) used in the following. 
Ax 4 + Bx 3 + Cx 2 + Dx + E = is a quartic if A ^ 0. A depressed quartic is a quartic such that 
B = holds. Ax 3 + Bx 2 + Cx + d — is a cubic equation if A ^ and in particular a cubic 
without x 2 term is called as depressed cubic. 

The first step in solving the quartic consists in transforming it into a depressed quartic. i.e. 
eliminate the x 3 term. Transform the equation into x 4 + ^x 3 + ^x 2 + ^x + M = and substitute 
x = u — jj. Simplifying the original quartic yields 

(5.1) u 4 + au 2 + Bu + j = 

where a = + §, / 3=^-^ + fand 7 = 5^ + T& ~ £p + f > and «1- O is a 

depressed quartic function, which is tantamount to 

(5.2) (u 2 + a) 2 + 0u + j = au 2 + a 2 

The next step is to insert a variable y into the perfect square on the left side of eq. (|5.2|) . Add 
both 

(u 2 + a + y) 2 - (u 2 + a) 2 = 2yu 2 + 2ya + y 2 
= (a + 2y)u 2 - 2yu 2 - au 2 

to ea. (|5.2p yields 

(5.3) (u 2 + a + y) 2 = (a + 2y)u 2 - (3u + (y 2 + 2ya + a 2 - 7) . 

The next step consists in substituting for y such that the right side of eq. (|5 .3[) becomes a square. 
Observe that (su + t) 2 = (s 2 )u 2 + (2st)u+ (t 2 ) holds for any s and t, and the relation between the 
coefficients of the rhs is (2st) 2 = A(s 2 )(t 2 ). Therefore to make the rhs of eq. (|5.3[) into a perfect 
square, the following equation must hold. 

B 2 

(5.4) 2y 3 + bay 2 + (4a 2 - 2-y)y + (a 3 - 07 - ^-) = • 

Similarly, transform eq (|5.4[) into a depressed cubic equation by substituting y = v — |a 

\ q , a 2 . , a 3 «7 B 2 . 

(5.5) w 3 + ( ___ 7)u+( __ + ^_P_ )=0 . 

Set P = — — 7 an d Q = — j^g + ^ — ■ Select any of the solutions of eq. (|5.5p of the form 



■§j — U where U = \ ^ ± -J ^- + S and {/ ^ 0, otherwise 1; = 0. In view of y = v — | this 
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solution yields for eq. (|5.4|) y = — |a + ^ — U for [/ ^ and y = — |q for [7 = 0. Now the rhs of 
eq. (|5.3[) becomes 



(a + 2y)?i 2 + (-/3)u + (y 2 + 2ya + a 2 - 7) = ^(ya + 2y) it + 
Combined with eq. (|5-3[) this yields two solutions for u: 



2^a + 2y 



±^/a + 2y± s j-[3a + 2y± 
u = 



2/3 



V<x+2y 



2 

where the first and the third ± must have the same sign. This allows us to obtain the solutions 
for x: 



ft 



±y/a + 2y±J- (3a + 2y± 



2/3 



B v a V V V^+2y 

X = -U + 2 • 

In particular for a; 4 - 5a; 3 - x 2 + 5x - 1 = 0, A = 1, B = -5, C = -1, D = 5, E = -1, and 
hence a = -f , (3 = -±^, P = -f and Q = |f§ and U w 1.21481 - 0.54955i 7^ 0, therefore 
2/ « 6.21621, and the solutions are p 3 ta 0.21982, C2 « 5.00829, C3 « -1.07392, C 4 « 0.84581. As for 
the equation a; 4 + 3a; 3 — a; 2 — 3a; — 1 = 0, the corresponding solutions are £5 ~ —0.53243 + 0.11951i, 
Ce ~ -0.53243 - 0.11951i, C7 ~ 1.10477 and Cs « -3.03992. 
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